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, To treat effects of electron correlations in geometrically frustrated pyrochlore and checkerboard 

fT^ ■ lattices, an extended single-orbital Hubbard model with nearest neighbor hopping ~ t and Coulomb 

repulsion ^-^ V is applied. Infinite on-site repulsion, U oo, is assumed, thus double occupancies 
of sites are forbidden completely in the present study. A variational Gutzwiller type approach is 
^ . extended to examine correlations due to short-range 1/— interaction and a cluster approximation is 

developed to evaluate a variational ground state energy of the system. Obtained analytically in a 
special case of quarter band filling appropriate to LiV2 04, the resulting simple expression describes 
^ ■ the ground state energy in the regime of intermediate and strong coupling V. Like in the Brinkman- 

Rice theory based on the standard Gutzwiller approach to the Hubbard model, the mean value of 
the kinetic energy is shown to be reduced strongly as the coupling V approaches a critical value Vc. 
This finding may contribute to explaining the observed heavy fermion behavior in LiV204. 
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I. INTRODUCTION 



The pyrochlore lattice is a three-dimensional network of corner-sharing tetrahedra formed by B-cations in py- 
rochlore oxides A2B2O7, in cubic spinel oxides AB2O4 and in cubic Laves-phase intermetaUic compounds AB2. These 
^ r compounds exhibit a wide variety of physical properties ranging from magnetic insulator through bad metal to su- 
00 ' perconductor with relatively high transition temperature. One of the most spectacular properties in this family of 
, d-systems is the heavy fermion (HF) behavior found in LiV2 04 and Y(Sc)Mn2 at low temperaturesiS*^ To under- 
stand the nature of HF quasiparticles in Y(Sc)Mn2, it was suggested'''^ that effects of geometrical frustration on the 
itinerant electron system in the pyrochlore lattice are of crucial importance. 

In fact, in the pyrochlore lattice with a single orbital on each site, a special lattice geometry manifests itself through 
(•~^ a flat band on the top of the dispersive one-particle electron spectrumj^ In the half-filled case, the Fermi level touches 
"Y^ ' the flat band singularity from below. In actual band structure calculations of a multi-orbital d-system, this singularity 
d is removed and, instead, a sharply peaked feature in the density of states is seen. By treating the problem of Y(Sc)Mn2 
, within a single band pyrochlore Hubbard model in a weak correlation limit and assuming proximity of the Fermi level 
I ' to the peaked feature in the density of states, the authorgii^ suggested a new scenario for the HF-like behavior of 
' ^ the frustrated system. In contrast, in LiV2 04 the mean occupancy of the V-ion d-shell equals 1.5 and corresponds 
to the quarter-filling of the almost three-fold degenerate t2g band.^ Therefore, in this compound the Fermi level is 
^ [ pushed down strongly from the 'flat-band' peculiarity and the above arguments for the HF quasiparticle formation 
• • ' are not applicable any more. In this respect, even by taking into account strong on-site Coulomb correlation effects, 
not included in the band structure calculations, one cannot cure the situation. 

Several scenarios for explaining HF properties in LiV2 04, including the one based on multi-component fluctuations^, 
^ ' are proposed in the literature and reviewed briefly in Ref0- While the strong on-site Coulomb interaction ~ [/ is 
5^ , incorporated necessarily in most of these scenarios, effects of inter-site Coulomb correlations on the HF properties in 
the pyrochlore lattice have been studied much less. The latter effects in the form of Anderson's 'tetrahedron rule- 
are used explicitly in ReflTol to build a physical picture of low temperature properties in LiV204. According to this 
picture. Coulomb repulsions of electrons on nearest neighbor sites suppress charge fluctuations appreciably and a state 
of slowly fluctuating and nearly decoupled chain-like and ring configurations of spins s = 1/2 is formed to give a large 
low temperature spin entropy. However, the mechanism for suppressing charge fluctuations in the pyrochlore lattice 
has not been elaborated in detail. 

In the present paper, we address this problem and concentrate on the study of charge degrees of freedom and their 
possible role in forming heavy quasiparticles. For this purpose, the potentially important role of inter-site electron 
Coulomb interaction in the pyrochlore lattice is emphasized in a particular case of quarter band filling appropriate to 
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LiV204- We use an extended single-orbital Hubbard model composed of the kinetic term with a tight-binding nearest 
neighbor (n.n.) hopping amplitude t, strong Hubbard term ~ U and Coulomb term describing n.n. electron repulsion 
^ V oi variable strength with respect to t. In spite of its serious simplification, this minimal model is believed to 
capture generic effects of strong inter-site Coulomb correlations inherent to more realistic pyrochlore electronic models 
including several d-orbitals on each lattice site. Noting that U is the largest energy parameter, we simply put U ~^ oo 
in our calculations. At t/ ^ cxd and far away from half-filling, one expects the system to be in a metallic state at 
weak coupling V. Our aim is to show that in the strong coupling limit, V ^ t, an effective electron bandwidth in 
the quarter-filled system can be strongly reduced due to short-range charge correlations. This regime of a strongly 
correlated metal can be understood as a state of almost localized electrons. 

Short-range spin correlations are neglected here because the limit ?7 — *■ oo is assumed. Instead, we concentrate 
completely on the problem of short-range charge correlations which are believed to provide the dominant mechanism 
for strong suppression of the electron kinetic motion in the pyrochlore lattice and, therefore, may contribute to the 
heavy quasiparticle formation. 

To conclude this Section, we outline briefly the structure of the paper and the method used. We extend the 
Gutzwiller variational approach to a form that allows us to treat the inter-site correlation problem. For completeness, 
we involve in our consideration also the 2D checkerboard lattice which is a planar analogue of the 3D pyrochlore lattice. 
Our method is applicable equally to both cases and relies on a cluster character of the pyrochlore (checkerboard) 
lattice structure, which means that each tetrahedron (plaquette) is regarded as an elementary entity. In Section II, 
the underlying electronic model is defined and rewritten in a cluster representation. The ground state variational 
wave function of the interacting system is constructed by applying to the system of noninteracting electrons two kinds 
of projectors. The first is the Gutzwiller projector taken in the limit t7 ^ oo to eliminate double occupancy of sites. 
The second projector is defined in terms of cluster occupancy operators in such a way that the weights of different 
charge configurations in each cluster are controlled by a correlation strength parameter rj which is the only variational 
parameter in our theory. 

In Section III, a variational ground state energy functional is defined and a factorization procedure of its approximate 
calculation is discussed in detail. Within this context, the concept of the basic cluster in the pyrochlore (checkerboard) 
lattice is introduced and a mechanism for suppressing of inter-cluster charge fiuctuations is explained. In Section IV, 
the calculations are performed and the resulting ground state energy as a function of the coupling V is analyzed. 
This analysis shows that the mean kinetic energy in the system is strongly renormalized and tends to zero as V 
approaches some critical value Vc- Also, our variational theory suggests, that electron localization transition occurs 
a.t V = Vc, and for V > Vc the electrons are in a fully localized insulating state for which quantum charge fiuctuations 
are completely missed. The results are discussed in Section V. By drawing a close analogy with the Brinkman-Rice 
localization transition predicted in the variational studj*ii of the standard Hubbard model, we argue that improper 
description of the actual metal-insulator transition does not invalidate our cluster variational theory in all respects. 
In particular, the theory is still valid in describing a metallic strongly correlated (the not yet localized) regime for 
V ^Vc- Concluding remarks can be also found in Section V. 



II. MODEL FORMULATION AND TRIAL WAVE FUNCTION 



Let us start with an extended single-orbital Hubbard model in the pyrochlore or checkerboard lattice 

H = -t [c\^Cja + c]„Cia) +"^^51 +^^71,^71,^, (1) 

<ij>,a <ij> i 

where the summation in the first two terms is taken over the n.n. sites, each pair < ij > is counted once; c|^(cja) is 
the creation (annihilation) electron operator with spin a =t,i and rii = ^„ c\^Cia- We are interested in the model 
regime far away from the band half-filling in the limit [/ — > oo. 

Let us consider a plaquette in the checkerboard or a tetrahedron in the pyrochlore lattices as an elementary /C-th 
cluster, each contains four sites numbered as fc = 1,...,4. Each lattice site belongs to two clusters and thus the 
neighboring clusters overlap. If L is the number of the lattice sites, then L/2 is the number of overlapping clusters, 
i.e. IC — 1,2,..., L/2. Now, the first two terms in the Hamiltonian JQ) can be rewritten in the cluster representation 

L/2 

Ht-v = ^ Hfc, Hjc = T/c + V/c, 

!C=1 

TK = -t J2 E4.cfe'., (2) 

k^k' <y 
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Since U oo, the double occupancy of sites is forbidden and the unit operator on each k-th site is 

Ik^ Pk{0) + Pk{l), 

Pfe(O) = (l-nfcT)(l-7ifcx), 

Pfc(l) = J2 ^fe-(l)' = nfc,(l - Uka), 



(3) 



where Pfc(O) and Pfe(l) are the projection operators onto the empty and singly occupied fc-th site states, respectively; 
a = —a. Further, to distinguish cluster states with occupancies n = 0, . . . , 4, we define the cluster projection operators 



VK{0) = \{Vk{0), 



k=l 
4 



fc=l 

P^(2)- ^P,(l)Pfe,(l) 

k<k' 
4 

fc=i 

4 

rK{^) = \{Fk{i), 



n P'^'(o) 



n 

k"^k,k' 



(4) 



n P'^'(i) 



k=l 



where {k,k',k") G /C is implied. Disregarding the spins, we note that the configurational space of a cluster contains 
in total 16 allowed charge configurations. For instance, the occupancy state n = 2 is composed of 6 configurations. 
Taking into account the frustration of the system, we suggest that the basic operators I0J and the associated cluster 
occupancies in the lattice suffice to capture the correlated metallic regime of the model, as suggested by the form © 
of Vjc term below. 

In terms of the projection operators the kinetic and the Coulomb terms in (j^J can be presented as follows 



(5) 



= V^[^yc(2) + 3P;c(3) + 6P;c(4)], (6) 

where the diagonal form of Tx, occurs because any electron hopping within a cluster does not change its charge 
occupancy. Finally, we define the cluster number operator TVjc — X^fceK;"-*:' note that the mean value (A^ac) — 2 
corresponds to the quarter band filling. 

A trial ground state wave function of the interacting system is constructed below in two steps. First, in the limit 
V — Q, we introduce the Gutzwiller projected Fermi sea [U — > oo): 

1 ^ 

\il^)v=o = ^ n [1 - "fcT^fcil l^o) = IV'o), (7) 
k=i 

where |$o) is a single Slater determinant describing a metallic state for N — L/2 noninteracting electrons; Vl is 
the norm. Since the system is away from half-filling, the Gutzwiller projectorsiSii^ in {Tj) are expected to lead to a 
moderate renormalization of metallic properties. A regime of strongly correlated metal at quarter band filling may 
occur in an interacting system with strong inter-site Coulomb repulsion V . 

To take into consideration electron correlations due to ^ ^ 0, we define below a product of generalized Gutzwiller 
projectors Gic{i]) acting in the occupancy space of the /C-th clusters (/C = 1, 2, . . . , L/2); here t] is the variational 
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parameter varying in the range < 77 < 1. To justify the proper choice of G)c{r]), let us start with the Hmit of extreme 
coupHng, V/t ^ 00. In this hmit, the ground state cluster charge configurations are those obeying the 'tetrahedron 
rule', i.e. each cluster is occupied strictly with two electrons, and the ground state energy is V per cluster. Therefore, 
Giciv) — ^k;(2) and our convention is that 77 = in this limit. At finite t, electron hopping processes create exited 
charge configurations of higher, n = 3,4 and lower, n = 0, 1, cluster occupancies. To ascribe optimal weights to the 
exited configurations, we use the following prescription. Consider first a pair of clusters, JCi and IC2, occupied with 
ni = 1 and 712 = 3 electrons. By noting that the potential energy of such a fluctuation is 3V, which exceeds the 
energy 2V of the dominant configuration of two electrons per cluster by V , we ascribe to it a weight factor ry < 1. 
The cost of the potential energy of the more extreme charge fluctuation of a pair of clusters, i.e. ni = and n2 = 4, is 
4:V, and therefore we associate the weight factor 77^ with this fluctuation. More complex many-cluster excited charge 
configurations can be easily recognized to be combinations of the elementary ones discussed above. These observations 
are comprised in the following form of a generalized cluster Gutzwiller projector 

Gk{v) = V^V^iO) + VvVKil) + Vk{2) + VvVd^) + V^V^i^), (8) 
and of a trial ground state wave function defined as 

L/2 

mv))^l[^K{vMo) = \i'), (9) 
/c 

where 0<77<1- If?7 = l7 the projector Gic{i]) becomes the unity operator, which corresponds to the limit of zero 
coupling, = 0. 

Note, the form of a generalized Gutzwiller projector and, hence, the trial wave function jnj are inspired by the 
form of Vic term and Anderson's reasoning^ on a macroscopical degeneracy of the ground state charge configurations 
in the limit of the extreme coupling, V/t — > 00. We expect the wave function ansatz @ to be valid down to intermediate 
coupling values V ^ zt, where 2; = 6 is the pyrochlore (checkerboard) lattice coordination number. It is much less 
appropriate in the limit of a weak coupling V t, where the kinetic energy effects become dominant and the arguments 
based on the potential energy counting only arc not sufficient. 



III. VARIATIONAL GROUND STATE ENERGY AND FACTORIZATION PROCEDURE 

With definitions and the ground state energy of the lattice system £{ri) — {ip\H\ip) / can be written 
as a sum of cluster energies 

Siv) - -,Ei,l Ei,) . ^^(^)|^(^)^ , (10) 

where the translational invariance of the wave function @, i.e. {ip\Hic\ip) = {ip\Hic'\'ip) for K, 7^ K.' , is taken into 
account. From now on, the expectation energy value calculation reduces to a local problem. 

To evaluate E{ri), an approximate procedure is developed and applied below . The starting point is to divide the 
lattice into two parts. The smaller one, which is refereed below to as the basic cluster (be), contains five connected 
elementary clusters, as explained in Fig.l for the checkerboard lattice. According to this picture, the same hopping 
parameter t and Coulomb interaction constant V should be ascribed to the horizontal, vertical and diagonal bonds 
(linear segments) connecting neighboring lattice sites. One can easily realize that the basic cluster in Fig.l has its 3D 
counterpart formed by five corner-sharing tetrahedra in the pyrochlore lattice, the former being a planar projection of 
the latter. Therefore, both cases can be treated on an equal footing. With this division, we suggest that the essential 
effects of the local interactions involved in E{ri) are captured by a basic cluster. 

Below, a factorization procedure for the trial wave function is used, which reduces the calculation of expectation 
value £{1]) from the lattice to a basic cluster. In the present study, we take advantage of another simplifying 
assumption which amounts to the neglect of the spatial charge correlations in the starting wave function Q. Such 
an approximation is inherent also to the Gutzwiller approach study of standard Hubbard model. Actually, as shown 
in Ref jlJi, when calculating the hopping matrix elements (the interacting U -term is accounted for exactly) in the A'^- 
particle state j^o) of noninteracting electrons, the neglect of the configuration dependence at this stage is equivalent 
to all assumptions made by Gutzwiller— and leads to an identical result.^™ The approximation is controlled by the 
parameter where z is the lattice coordination number; it gives an exact result in the limit of infinite spatial 
dimension and requires only very small corrections for three-dimensional lattices^^ii^ Though the wave function Q 
describes an electron system with extremely strong on-site repulsion, we expect that foregoing arguments are still 
valid provided the electron concentration is sufficiently far away from the half filled band case. 
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After explaining the very essence of the approximation chosen to calculate the expectation energy value E{ri), 
we proceed now with its evaluation. First, a basic cluster is selected as a lattice fragment composed of a central 
elementary 0-cluster whose four sites are shared with four (I,...,IV) side elementary clusters, as depicted in Fig.l. 
From now on, the trial state @ can be written in the following form 

mv)) = Gko(^)G{k}(77)G{;c'}(^)IV'o), (11) 

where Gtc^ is referred to the central 0-cluster, Gsjc] ~ 11x^=7 is for four side clusters and Gj/c'} — ^k'{-^Kq k) G/C' 
is for the remaining system. The factorization procedure now reads as follows 

(7Ao|G'{K:}GK:ogK:oGK:oG{K:}|V'o)(V'o|G^^,}|V^o) _ 
^"^^ ^ (V'o|G2^jG2^^j|^o)(^o|G2^,j|V'o> 

_ (V'o|G{K:}GK;Qi?/CoG/CoG{/c}lV'o) (V'&cI^^KoIV'&c) 



(12) 



(V'olG^^jG^^^jlV'o) (V'bclV'bc) 
where the trial state \t\}hc) in the cluster approximation is taken in the form 

|V'6c(?7)) =GKo(r/)G{K}(77)|V'o), (13) 

and the labelling he is due to a distinguished basic cluster either in the pyrochlore or in the checkerboard lattice. 
Both, three- and two-dimensional basic clusters have the same lattice connectivity and, therefore, in our study 
are treated equally. We believe that the basic cluster chosen is inevitably a minimal one for the present problem. 
This means that it could not be reduced further, for instance, to the size of 4-site elementary cluster, otherwise 
the most important correlation effect of the kinetic energy reduction would be lost. Actually, when calculating 
('0bc|I)c:ol''/'6c)/('0bc|'0bc)/ = {TKf^bc, One finds that the electron hopping processes involved in T^Co do not change 
an occupancy of the /Co -th cluster, but any of them changes simultaneously occupation numbers of a pair of side 
clusters. As a result, the side cluster charges deviate from their mean value n = 2. With increasing coupling constant 
however, the cluster states of high, n = 3,4 and low, n = 0, 1, occupancies are suppressed, and thus the hopping 
processes within the JC^-ih cluster are hampered. As a consequence, the expectation value of the kinetic term {T]Cg)bc 
is reduced strongly in a way similar to that found in the standard Hubbard model treated within the Gutzwiller 
approximation. In the latter case, a minimal basic cluster is 2-site one as suggested by Razafimandimby.^^ 

In order to facilitate evaluation of E(ri), we resort to a new diagonal form of the operator product Gic„G^^i<^j 
with respect to the occupancy number n of the central cluster ICq. For this purpose, let us consider side clusters 
K. — I, . . . ,IV and their labels K. as composite ones, /C = (fc, k), for instance, I = (1, 1), etc. Here k denotes a site 
common both to the /C-th and the central /Co clusters, while k refers to the complementary 3-site fragment of the /C— th 
side cluster, as shown in Fig. 1. For the fc— th fragment, let \k,m) be a quantum state with an electron occupancy 
m{= 0, .., 3) and V-^im) = \k, m) {k, m\ are the corresponding projection operators. In terms of Vjim), the projection 
operators Vic{n) of side /C-clusters can be now written in a split form V/cin) = Pk{0)'P-i:{n') + Pfc(l)Pj:(n — 1), with 
the convention 7'j7(— 1) = 'Pj:(4) = assumed. With this redefinition of V/cin), the projectors G/c of side clusters are 
also split in two terms, each selecting one of two possible occupancies, n = 0, 1, of the fc— th site common to a given 
/C and the central /Co clusters: 

Gk(.v) = PmQuiv) + Pk{l)R^{il)- (14) 
Here Q-^ and R-^ operate in the charge configuration space of the fc— th fragment 

Rliv) - V^7'^(0)+P^(l) + V^7'^(2)+ ,72^^(3). (15) 

with the following obvious properties, Q^iv) = Qk^V^) ^"^"^ ^^^^ ~ ^k^^^^- 

It is not a difficult task to rewrite the operators V-j:{m) and their linear combinations H15|l through the site projection 
operators ||3J), which is useful for further calculations. The weight factors in (|15|l serve as a measure of the correlation 
strength in dependence on the fc— th fragment charge occupancy. Finally, we use the definition © for GjCj, and the 
presentation lfTl|l for G^; to write down the following form of the product GjCgG^ic}: 

4 

GiCo{v)G{K}{v) = 5]/„(ry)X(n;r;), (16) 

n=0 
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where weight factors fn{il) a-re found to be fo = fi = rf, /i = /s = y/Vy /2 = 1, and the operators X{n; rf) discriminate 
the basic chister charge configurations in accord with the central cluster occupancy n: 

4 

X(0;77) = X{PmQki^), 

k=l 

4 4 

X{l-rj) - ^Pfe(l)%(,7) II P,,{0)Q-{rj), (17) 

k=l k'^k 

X(2;77) - ^Pfe(l)Pfe,(l)%(r,)i?^(ry) [] P^„ {Q)Q-{^r^) ^ 

k<k' k"iLk,k' 
4 4 

fc=i fe'j^fc 

4 

X{A-ri) = n^fc (!)%('?)• 

fe=i 

The factorization procedure (|12|l constitutes the first step of our approximation for calculating the expectation en- 
ergy value E{ri). As mentioned before, the second approximation consists in the neglect of electron spatial correlations 
within the Gutzwiller projected Fermi sea state jV'o) — \4')v=0: Eq.jTI). For instance, when having a two-site electron 
density correlation function (■(/'ol^-fc.o-'T-fc',^' |'0o)i we decouple it as follows (fc ^ k'): 

(V'o|"-fc,<T»^fc',<T'|'0o) - {^Jo\nk,a\'^o) {^Jo\nk\a'\i'o) = Q") (^^) 

where n = N/L is an electron concentration and the extra factor 1/2 in the right-hand side is due to spatially 
uncorrelated electron spins. Since doubly occupied sites are forbidden, ('0o|^fe,T^fe,il'0o) = 0; in terms of site projection 
operators (|3J), the decouphng lfTH|) reads (fc ^ fc'): 

(V'o|Pfe(^)^^c'("')l^o> - {Pkin))o{Pk'in'))o, (19) 

where {Pk{n))o = (^/'ol-Ffc('^)lV'o) = l/2('5,i.o + <5n,i)- Multi-site density correlation functions are decoupled in the same 
manner. 

To complete this Section, let us calculate the norm 

4 

where the average (...)o is over the state \ipo). 

By using the approximation (|18|l and p9(l for a given electron concentration n =^ 1/2, one finds an intermediate 
result 

{Xin;v'))o = (^ly C^A{Q{v'))o)'~'' mv'))o)\ (21) 

where is a binomial and both {Q-^}o and {R-^)o are calculated to be fc— independent expectation values {Q{'rj^))o = 
(i?(r?2))o = (^4 + 47y + 3)/8. 

From Eq.lEOI), the norm can be now written as 



(V'f)c|V'bc) 



77^ + 4?7 + 3' 
8 



(22) 



IV. CALCULATION OF THE VARIATIONAL ENERGY 



We start with the evaluation of an expectation value of the Coulomb energy. Since the procedure is similar to that 
used in the previous Section, we present the resulting expression for the mean interaction energy per cluster 

{^bc\VKo\i'bc) ^ ^yt+^h+l ^ H^' ^ = ^ .23) 

(V'iclV'bc) 774 + 477 + 3 I ,^ = 
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Both limiting values in (|23|l can be easily understood on a general ground. Indeed, the value {3/2)V corresponds to 
spatially uncorrelated electrons homogeneously distributed over the lattice, while the value V is in compliance with 
any lattice charge configuration having two electrons within each elementary cluster. 

The evaluation of the kinetic energy is much more involved. As a starting point, one has to rely on a general 
expression 

k^k' (J 

3 

= E E E yL^k-aX{n- ))o, (24) 

k^k' a n=l 

where each pair of sites (fc, k') G K-o contributes equally to the sum l|24() : from now on we also use shortened notation 
X{n;ri) = X{n;). To make the calculation more transparent, we introduce projected fermionic operators cj^^(cfc(j) 
related to cl^icka) as follows 

4. = Pk{l)clPkiO), Cka = Pk{0)ckaPk{l). (25) 

When operating in a space with no double site occupancy, cl,^{cka) and c[^(cfe(j) are equivalent. 

Let us consider a particular pair (fc = 1, A:' = 2) and insert into H24|l the explicit form of X{n; ) coming from H17|l . 
Then we obtain successively for n = 1,2,3 the following intermediate form of hopping matrix elements 

= {R-iclQi ■ Q-2haR-2 ■ PmQl ■ A(0)Ql)o, 

(X(2;)cLc2aX(2;))o= (26) 

- {R-AaQl-Q2i2aR-2 ' [^^3 (0)Qi • P4 (l)i?l + ^3 (l)i?i ' ^4 (0)Q1] )o , 

(X(3;)4c2.X(3;))o = 

= {Rlc\aQi ■ Q-2C2.R2 ■ P3{l)Rl ■ PiimDa- 

For brevity, in the above expressions we drop the arguments in Qj^iff) and _Rj.(r/). 

It is worth discussing shortly the physical content of expressions (I26II . An electron hopping within the central 
cluster is affected by its charge surrounding, which may occur in many configurations. In our description, all allowed 
configurations of the charge surrounding are comprised by a set of operators Qj. and i?j involved in the hopping 
matrix elements (|^ . This becomes more clear if one uses the explicit form of Qj. and Rj., Ea. ((T5|l . which breaks up 
a hopping amplitude from H2t)|) into a sum of terms, each of them corresponds to a particular configuration weighted 
with some 77— dependent factor. In the uncorrelated limit, 77 = 1, all these factors are unity, while for 77 < 1 hopping 
amplitudes are suppressed due to short-range charge correlations. 

The neglect of spatial electron correlations in |'!/;)o allows us to approximate the amplitudes H26|l by using a decou- 
pling procedure. For instance, the first expression in is decoupled as follows 

(X(l;)cLc2.X(l;))oC^ 

~ {RiQi)o{R-2Q-2)o{Ql)o{Q^MPM)o{Pm)o{A.~^2.)o = 

- (1/2)2 (i?Q>2(g2)2(ct^C2.)o, (27) 
and the remaining two amplitudes in H26|) are approximated in the same way. In (|27|) the last equality takes into account 
fc— independence of expectation values {Ri.Q^)q = {RQ)(), {Q^o = {Q'^)o and the property {c[^C2a)o = {cia-C2a)o- 
The expression for {Q^ {v))o = {Q {v'^))o, {R^ {v})o — {R {v^))o, is given in the previous Section and the expectation 
value {Rf. (r/) Qj. {ri))o — {RQ)o is calculated easily to give {RQ)o = {v^ + 3) /4. 

Next, by collecting these results and with the reference to H24() . we obtain the following relation between hopping 
amplitudes in the correlated, V 0, and 'uncorrelated', V — 0, states 

(cLcfcv)ftc ^^Viv + 1) iv' + 3)' iv' + + 3)' (cLcfcv)o, (28) 

The latter hopping amplitude {cl,^Ck'a)o is obviously related to the expectation value of the kinetic energy per cluster 
calculated at V = 0: 

{TKo)o^-t (4.Cfc'.)o. (29) 

{k,k')^Ko 
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This leads to the final result 

mKM)bc . , , , 167? (7; +1) (77^ + 3)^ 

wmbc (774 + 47; + s)"* 

where q{ri) is the kinetic energy renormalization factor. Note that in the 'uncorrelated' limit, V = 0, the necessary 
condition (7(7/ = 1) = 1 is fulfilled. The kinetic energy H3U|) and the Coulomb term H23|) together give the final result 
for the variational energy E{ri) per cluster. 

In the extreme limit, V/t 00, all clusters are doubly occupied (the 'tetrahedron rule' is fulfilled completely), 
which leads to the upper bound of the variational energy, Emax — V , at i] = 0. Let us consider 



(774 + 477 + 3) 



V-8\£\ 



(772 -77 + 1) (774 + 477 + 3)' 



< 0, (31) 



where {T]Cg)a = —\£\ < 0. The expression between the brackets in (|31|) changes sign at some critical coupling V — Vc, 
with the following requirement 

A close inspection of Eq. (|31|l suggests the critical value Vc — 8\£\. To check this suggestion, let us minimize £^(77) by 
solving the equation dE{rj)/drj — 0, which relates 77 to y and requires the only solution r] — atV — Vc — 8|£|. In 
the vicinity of V^ i.e. for V < 8|£|, the kinetic energy renormalization factor is given by 

which is checked numerically to describe the mean kinetic energy reduction with high accuracy in a wide range of 
coupling, 1/2 < V/Vc < 1. Rather generally, one may estimate the value of {Tic„)o to be £ = —azt, where the 
coefficient a is of order unity. Therefore, according to our convention, the solution H33I) applies to the strong coupling 
regime of the model (^1, where the wave function ansatz in the form @ is valid. We do not intend to discuss here 
the variational ground state energy in the weak coupling limit, V ^ t. As argued in Sections II, in this limit kinetic 
energy effects dominate inter-site Coulomb correlations, the wave function ansatz in the form jnj is not applicable 
any more and its complementary variational extension is required. 

The physical picture emerging from these results is discussed in the following Section. The discussion is closely 
related to results of the Brinkman-Rice theory>ii 



V. DISCUSSION. 



Based on a simple trial wave function (|5J) with a single variational parameter 77, we infer from the results (|30|l . H31|) 
and H33|) that in the strong coupling regime, V ^ t, the mean kinetic energy per cluster, 9(21co)oj of the quarter-filled 
pyrochlore electronic system is strongly reduced by the factor q {—^ 0, as — > — 0). This reduction means the 
increasing electron localization in real space and thus a strongly correlated metallic state is expected to appear for 
V ;$ Vc- For V > Vc, the mean kinetic energy is zero and the variational energy is V per cluster, which describes 
an insulating highly degenerate state of fully localized electrons. Such a description of the insulating state is correct, 
however, only in the limit of extreme coupling, V/t — > cxd. For finite, however, still large V/t, both in the pyrochlore 
and checkerboard latticesii charge fluctuations to leading order in t/V would lower the ground state energy by the 
value of order t'^/V. It means that even in an insulating state, full electron localization cannot be reached. This 
observation casts doubt on the very existence of the predicted localization transition. 

In this respect, our results are very similar to those obtained by Brinkman and Rice in their stud}>ii of the Hubbard 
model based on Gutzwiller's wave function and approximation in the limit of a half-filled band. There, the driving 
force for electron correlations and localization, at a finite strength Uc, is the strong on-site Hubbard interaction U, 
while in the present study of the quarter band filling the electron localization is due to inter-site Coulomb interaction 
V. The merits and failures of Brinkman-Rice theory are well understood and widely discussed in the literature (see 
Refs. [1I1IIT9II20L and references therein). It was realized, for instance, that the simple Gutzwiller wave functiorki* 
is not rich enough to describe the true insulating state in the Hubbard model and that at any finite dimension of a 
lattice the occurrence of the Brinkman-Rice localization transition is merely due to the Gutzwiller approximation used. 
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According to this approximation, the spatial correlations in the system of noninteracting electrons are neglected. (In 
the present study we neglect spatial correlations in the Gutzwiller projected Fermi sea Q which is the starting wave 
function chosen before applying a generalized Gutzwiller projector (jHJ). This finding, however, does not invalidate 
the general conclusion about the usefulness of the Gutzwiller approximation which is known to give a proper physical 
description in a number of situations. Indeed, as shown by VoUhardt, Wolfle and Andcrson^-'^, a Hubbard lattice-gas 
model in the not yet localized regime {U/Uc ^ 1) describes ground state properties of normal liquid ■^He ('almost 
localized' Fermi liquid^^). 

In view of the above discussion, our results suggest that in the quarter-filled pyrochlore lattice a metallic state of 
nearly localized electrons may occur mainly due to strong enough intersite Coulomb interaction, i.e. for V/Vc ^ 1. 
To our knowledge, the present variational theory is the first one showing the importance of short-range Coulomb 
interaction in forming a strongly correlated metallic state in the pyrochlore lattice. Calculation of physical properties 
of such a state is, however, beyond the scope of the present study. Therefore, many questions remain to be answered. 
For instance, whether the predicted state is the Fermi liquid and, if it is so, may one relate the renormalization 
parameter q to the discontinuity at the Fermi vector of the momentum distribution function? An answer would allow 
us to connect the inverse of g to a charge quasiparticle mass enhancement in the strongly correlated regime of the 
model. Concerning a description of the weak interaction case, V ^ t, the trial wave function (|5J| must be extended, for 
instance, by introducing a second variational parameter, which would make the present theory more attractive. We 
believe, however, such an extension of does not change our main results, and retain these problems for a further 
study. 
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FIG. 1: Fragment of the checkerboard lattice chosen to be a minimal basic cluster in the variational cluster calculations. Each 

label K, of side clusters (I, . . . , IV) is a composite one, K. — (fc, fc), where fc {= 1, 4) is referred to a site shared with the central 
0-cluster and k {— 1, 4) is for the complementary 3-site fragment (hatched). In the pyrochlore lattice, the same labelling is 
used for a minimal basic cluster formed by five corner-sharing tetrahedra. 



